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Abstract 

The SLIM experiment was a large array of nuclear track detectors located at the Chacaltaya high 
altitude Laboratory (5230 m a.s.L). The detector was in particular sensitive to Intermediate Mass Mag- 
netic Monopoles, with masses 10 5 < Mm < 10 12 GeV. From the analysis of the full detector exposed for 
more than 4 years a flux upper limit of 1.3 • 10" 15 cm~ 2 s" 1 sr~ 4 for downgoing fast Intermediate Mass 
Monopoles was established at the 90% C.L. 

1 Introduction 

In 1931 Dirac introduced Magnetic Monopoles (MMs) in order to explain the quantization of the electric 
charge, obtaining the formula eg = nhc/2, from which g = ngjj = nhc/2e = n 68. 5e = n 3.29 • 10 -8 in 
the c.g.s. symmetric system of units pQ; n is an integer, n = 1, 2, 3, ... MMs possessing an electric charge 
and bound systems of a magnetic monopole with an atomic nucleus are called dyons. An extensive 
bibliography on MMs is given in ref. [2]. Relatively low mass classical Dirac monopoles have been 
searched for at high energy accelerators 0] . 

Magnetic Monopoles are present in a variety of unified gauge models with a wide range of masses. 

Grand Unified Theories (GUT) of the strong and electroweak interactions at the mass scale Mq ~ 
10 14 -7- 10 15 GeV predict the existence of magnetic monopoles, produced in the early Universe at the end 
of the GUT epoch, with very large masses, Mm > 10 16 GeV. Such monopoles cannot be produced with 
existing accelerators, nor with any foreseen for the future. In the past, GUT poles were searched for in 
the cosmic radiation. These poles are characterized by low velocities and relatively large energy losses 
[5]. The MACRO experiment set the best limits on GUT MMs with g — gD, 2gD, 3gn and dyons at the 
level of - 1.4 • 10~ 16 cm~ 2 s" 1 sr" 1 for 4 • 10~ 5 < (3 = v/c < 0.7 0. 

Some GUT models and some supersymmetric models predict Intermediate Mass Monopoles (IMMs) 
with masses 10 5 < Mm < 10 12 GeV and with magnetic charges of multiples of go', these MMs may 
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have been produced in later phase transitions in the early Universe and could be present in the cosmic 
radiation [2 IB]. 

IMMs may be relativistic since they could be accelerated to high velocities in one coherent domain of 
the galactic magnetic field. In this case one would have to look for downgoing, fast (/? > 0.03), heavily 
ionizing MMsQ. 

The main purpose of the SLIM (Search for Light Monopoles) experiment at the Chacaltaya laboratory 
in Bolivia at 5230 m a.s.L, was the search for IMMs \9\. An exposure at a high altitude laboratory allows 
to search for MMs of lower masses, higher magnetic charges and lower velocities, see Fig. [1] 

The searches for IMMs by Earth based detectors are essentially limited to downgoing particles [TU] , 
Water Cherenkov detectors are limited to fast downgoing IMMs (with j3 > 0.5), and a search can be done 
if the detectors are able to discriminate against the large background of cosmic ray muons [llj . 




Figure 1: Accessible regions (above lines) in the plane (mass, f3) for monopoles with magnetic charge 
9 — 9d coming from above for an experiment at altitudes of 20000 m, 5230 m, and for an underground 
detector at the Gran Sasso Lab. (average rock overburden of 3700 m.w.e.) 



The SLIM detector was also sensitive to Strange Quark Matter nuggets [HI H3] and Q-balls [33] . The 
results on these Dark Matter candidates are discussed in ref. [T5] . 

In the following, we present a short description of the SLIM apparatus, the calibrations of the Nuclear 
Track Detectors (NTDs), the etching and analysis procedures, and the limits obtained by the experiment 
on IMMs and GUT Magnetic Monopoles. 

1 The interest in MMs was also connected with the possibility that they could yield the highest energy cosmic rays [8]. 
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Figure 2: Left: the SLIM modules installed at Chacaltaya. Right: composition of one of the 7410 
modules; each module was enclosed in an aluminized mylar bag filled with dry air at a pressure of 1 
bar. 



2 Experimental procedure 

The SLIM experiment was an array of NTDs d with a total surface area slightly greater than 400 m 2 [3] . 
The array was organized into 7410 modules, each of area 24 x 24 cm 2 . All modules were made up of: 
three layers of CR39® |, each 1.4 mm thick; 3 layers of Makrofol DE® 0, each 0.48 mm thick; 2 layers of 
Lexan each 0.25 mm thick and one layer of aluminum absorber 1 mm thick (see Fig. [2] right). The CR39 
used in about 90% of the modules (377 m 2 ) was of the same type used in the MACRO experiment [6]. 
The remaining modules, 50 m 2 , utilized CR39 containing 0.1% of DOP additive, CR39(DOP). 

Each module (stack) was sealed in an aluminized plastic bag (125 /im thick) filled with dry air at a 
pressure of 1 bar. The modules were transported to La Paz, Bolivia, from Italy in wooden boxes and 
their position with respect to the other modules in the shipping crate was recorded. The stacks were 
deployed under the roof of the Chacaltaya Laboratory, roughly 4 m above ground (see Fig. [2] left). The 
installation of the SLIM detectors started in February 2000 and ended in February 2002. The return of 
the material to Italy was organized in batches, after the completion of the 4 years exposure. 

The atmospheric pressure at Chacaltaya is about 0.5 bar; before shipping to Chacaltaya, in Bologna 
we checked the air tightness of the envelopes sealed with air at a pressure of 1 bar by placing a sample of 
them in an airtight tank at a pressure of 0.3 atm for a few months; no significant leakage was detected. 

From the experience gained with the MACRO Nuclear Track Subdetector [B] , we know that the used 
CR39 does not suffer from "aging" or "fading" effects for exposure times as long as 10 years [TB]. Further 
calibrations with 1 AGeV Fe 26+ ions in 1999 and 2005 and with 158 AGeV In 49+ in 2003 confirmed the 
quality and the stability of the CR39 used in the SLIM experiment [T7] . 

2 Another 100 m 2 of NTDs were installed at Koksil (Pakistan, 4275 m a.s.l.) since 2002 and were not used in the present 
analysis. 

3 The SLIM CR39 was produced by the Intercast Europe Co, Parma, Italy according to our specifications. 
4 Manufacturcd by Bayer AG, Leverkusen, Germany. 
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2.1 Environmental measurements 

During the first phases of the detector deployment we evaluated possible effects of climatic conditions on 
the detector response and possible backgrounds. Previous tests had shown that the CR39 response does 
not depend on the time elapsed from its production and the passage of the particle if the ambient tem- 
perature ranges between -20 °C and +30 °C. The minimum and maximum values of the air temperature 
in each detector hall in Chacaltaya was recorded 3 times a day over the lifetime of the experiment. The 
temperature values usually ranged from °C to 30 °C with an average value of 12 °C for the whole year 
and from one year to the other; however in the summer months in very few cases temperatures down 
to -5 °C were measured in the early morning. Therefore, no significant variations were expected in the 
detector response over the exposure period. 

We performed measurements of the radon concentration in different locations of the experimental 
rooms where the SLIM detectors were placed. We used for this purpose E-PERM® radon dosimeters. 
The measured radon activity was about 40 -j- 50 Bq/m 3 of air. According to our previous experience 
with the MACRO NTDs, we concluded that this level of radon induced radioactivity did not present a 
problem for the experiment, even in case of radon diffusion into the module bags. 

Two different types of neutron detectors (BTI bubble counters and a BF3 counter detectors) were used 
to measure the neutron flux at Chacaltaya, during the first installation shift of 2001 over the energy range 
of a few hundred keV to about 20 MeV [TS] ■ Neutrons of these energies interacting inside the detectors 
could induce background tracks, and their density could affect the scanning speed and efficiency. Both 
types of neutron detectors measured the accumulated dose. Consistent results were obtained by both 
types of detectors. The accumulated dose measured in open air and near the detectors was very similar. 
The absolute neutron flux was computed using the BTI bubble counters for which the efficiency is known. 
A value of (1.7 ± 0.8) • 10 -2 cm~ 2 s^ 1 was obtained, which is in agreement with other reported neutron 
flux data at the altitude of Chacaltaya and with more recent measurements at the same location |19j . 
The necessity to reduce the neutron induced background in CR39 required us to study special etching 
procedures, mainly based on the addition of ethyl alcohol to the etching solutions. As discussed in the 
next section, the addition of alcohol reduces the background tracks on the detector sheets and improves 
the surface quality (i.e. greater transparency), at the expense of a higher threshold [17] . 

2.2 Etching procedures 

The passage of a magnetic monopole in NTDs, such as CR39, is expected to cause structural line damage 
in the polymer (forming the so called "latent track"). Since IMMs have a constant energy loss through 
the stacks, the subsequent chemical etching should result in collinear etch-pit cones of equal size on both 
faces of each detector sheet. In order to increase the detector "signal to noise" ratio different etching 
conditions [El [17] were defined. The so-called "strong etching" technique allows better surface quality 
and larger post-etched cones to be obtained. This makes etch pits easier to detect under visual scanning. 
Strong etching was used to analyze the top-most CR39 sheet in each module. "Soft etching" was applied 
to the other CR39 layers in a module if a candidate track was found after the first scan. This process 
allows to proceed in several etching steps and study the formation of the post-etched cones. 

For CR39 and CR39(DOP) the strong etching conditions were: 8N KOH + 1.5% ethyl alcohol at 
75 °C for 30 hours. The bulk etching velocities were vb = 7.2 ± 0.4 /xm/h and vb = 5.9 ± 0.3 /xm/h for 
CR39 and CR39(DOP), respectively. 

The soft etching conditions were 6N NaOH + 1% ethyl alcohol at 70 °C for 40 hours for CR39 and 
CR39(DOP). The bulk etching rates were v B = 1.25 ± 0.02 /xm/h and v B = 0.98 ± 0.02 /xm/h for CR39 
and CR39(DOP), respectively. 

Makrofol NTDs were etched in 6N KOH + 20% ethyl alcohol at 50 °C for 10 hours; the bulk etch 
velocity was vb = 3.4 /xm/h. 
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Figure 3: Calibrations of CR39 nuclear track detectors with 158 A GeV In 49+ ions and their nuclear 
fragments with decreasing charge. The base areas (1 pixel 2 = 0.3 fim 2 ) of the etched cones were 
averages over 2 faces. The CR39 was etched in (a) soft and (b) strong etching conditions. 



2.3 NTD calibrations 

The CR39 and Makrofol nuclear track detectors were calibrated with 158 A GeV In 49+ and Pb 82+ beams 
at the CERN SPS and 1 A GeV Fe 26+ at the Brookhaven National Laboratory (BNL) Alternating 
Gradient Synchrotron (AGS). The calibration layout was a standard one with a fragmentation target 
and CR39 (plus Makrofol) NTDs in front of and behind the target [501 . The detector sheets behind the 
target detected both primary ions and nuclear fragments of decreasing charge. 

We recall that the formation of etch-pit cones ("tracks") in NTDs is regulated by the bulk etching rate, 
vb, and the track etching rate, vt, i.e. the velocities at which the undamaged and damaged materials 
(along the particle trajectory), are etched out. Etch-pit cones are formed if vt > Vb- The response of 
the CR39 detector is measured by the etching rate ratio p = vt /vb ■ 

After etching the standard calibration procedure was the following: 

(?) measure the base area of each track in NTDs with an automatic image analyzer system [21]. The 
projectile fragments carry the same (3 and approximately the same direction of the incident ion; the Z of 
each resolved peak is identified via the base area spectrum. The average base area distributions of the 
In 49+ ions and of their fragments in CR39, etched in soft or strong conditions, are shown in Figs. [3K,b 
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Figure 4: Reduced track etch rate (p - 1) vs REL for the CR39 (left) and CR39(DOP) (right) 
detectors, exposed to the 158 A GeV indium ion beam, etched in soft and strong etching conditions. 



(1 pixel 2 = 0.3 fim 2 ). 

(ii) For each calibration peak the Z/f3 is obtained and the reduced etch rate (p — 1) is computed. The 
Restricted Energy Loss (REL) due to ionization and nuclear scattering is evaluated, thus arriving to the 
calibration data of (p — 1) vs REL shown in Fig. 0]for both strong and soft etching conditions for CR39 
and CR39(DOP). For soft etching the threshold in CR39 is at Z//3 ~ 7 corresponding to REL ~ 50 MeV 
cm 2 g _1 . For strong etching the threshold is at Z/j3 <~ 14, corresponding to REL ~ 200 MeV cm 2 g _1 . 
The extrapolation of the calibration curves to p — 1 gives REL < 40 MeV cm 2 g _1 for soft etching and 
REL < 160 MeV cm 2 g _1 for strong etching. For CR39(DOP) the threshold in soft etching conditions is 
at Z/ (3 ~ 13 corresponding to REL ~ 170 MeV cm 2 g _1 ; the threshold in strong etching conditions is 
at Zjj3 ~ 21 corresponding to REL ~ 460 MeV cm 2 g _1 . The extrapolation of the calibration curves to 
p = 1 gives REL < 240 MeV cm 2 g _1 for strong etching. 

For magnetic monopoles with g = gjj, 2g£>, 3gn we computed the REL as a function of (3 taking into 
account electronic and nuclear energy losses, see Fig. [S] [2"2"] . 

With the used etching conditions, the CR39 allows the detection of (i) MMs with g = go for (3 ~ 10~ 4 
and for (3 > 10~ 2 ; (ii) MMs with g — 2gr> for (3 around 10 -4 and for (3 > 4 • 10~ 3 ; (Hi) the whole /3-range 
of4-10~ 5 </3<lis accessible for MMs with g > 2go and for dyons. 

For the Makrofol polycarbonate the detection threshold is at Z/j3 ~ 50 and REL ~ 2.5 GeV cm 2 g _1 
[TT] : for this reason the use of Makrofol is restricted to the search for fast MMs. 

2.4 Analysis 

After exposure at Chacaltaya the modules were brought back by air flights to Italy in order to be etched 
and analyzed in the Bologna laboratory. Three "reference" holes of 2 mm diameter were drilled in each 
module with a precision machine (the hole locations were defined to within 100 /zm). This allowed us to 
follow the passage of a "candidate" through the stack. The bags (envelopes) were opened, the detectors 
were labeled and their thicknesses were measured, using a micrometer, in 9 uniformly distributed points 
on the foil surface. 

The analysis of a SLIM module started by etching the uppermost CR39 sheet using strong conditions 
in order to reduce the CR39 thickness from 1.4 mm to ~ 0.9 mm. After the strong etching, the CR39 
sheet was scanned twice, with a stereo microscope, by different operators, with a 3x magnification optical 
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Figure 5: REL vs beta for magnetic monopoles with g — go, 2gD, 3gD- The dashed lines represent 
the CR39 thresholds in soft and strong etching conditions and the Makrofol threshold (see Sect. 12, 2 [I . 



lens, looking for any possible correspondence of etch pits on the two opposite surfaces. The measured 
single scan efficiency was about 99%; thus the double scan guarantees an efficiency of ~ 100% for finding 
a possible signal. 

Further observation of a "suspicious correspondence" was made with an optical 20 -j- 40 x stereo 
microscope and classified either as a defect or a candidate track. This latter was then examined by an 
optical microscope with 6.3 Q fc x 25 oc magnification and the axes of the base-cone ellipses in the front and 
back sides were measured. 

A track was defined as a "candidate" if the computed p and incident angle 6 on the front and back 
sides were equal to within 20%. For each candidate the azimuth angle tp and its position P referred to 
the fiducial marks were also determined. The uncertainties A9, Atp and AP defined a "coincidence" area 
(< 0.5 cm 2 ) around the candidate expected position in the other layers, as shown in Fig. 

In this case the lowermost CR39 layer was etched in soft etching conditions, and an accurate scan 
under an optical microscope with high magnification (500 x or 1000 x) was performed in a square region 
around the candidate expected position, which included the "coincidence" area. If a two-fold coincidence 
was detected, the CR39 middle layer was also analyzed. 

The bottom CR39 sheet was etched in about 50 cases; the third CR39 sheet was etched only in few 
cases, when there was still a possible uncertainty, and for checks (~ 16 times). Some Makrofol foils were 
etched for reasons similar to the previous point and for other checks concerning the Makrofol itself (~ 12 
times). 
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Figure 6: Illustration of the procedure used to define the "confidence" area where the possible 
continuation of a candidate track inside two (or more) sheets of the same module was searched for 
(see text for details). 



3 Results 

From the detector calibration we computed the SLIM acceptance for downgoing IMMs with g — gu, 2gn, 3gn 
and for dyons. For the i th module of area Si the acceptance was computed as 



V P / 

The total acceptance is the sum of all the individual contributions. 

Since no candidates were found, the 90% C.L. upper limit for a downgoing flux of IMMs and for dyons 
was computed as 



where At is the mean exposure time (4.22 y), SQ is the total acceptance, e is the scanning efficiency 
estimated to be ~ 1. 

The global 90% C.L. upper limits for the flux of downgoing IMMs and dyons with velocities j3 > 4- 10 -5 
were computed, as shown in Fig. [7J The flux limit for (3 > 0.03 is ~ 1.3 • 10 -15 cm~ 2 s _1 sr _1 . 

Two "strange events" were observed and were finally classified as manufacturing defects in a small 
subset of CR39 NTDs. These "strange events" are discussed in detail elsewhere [23] . 

4 Conclusions 

We etched and analyzed 427 m 2 of CR39, with an average exposure time of 4.22 years. No candidate 
passed the search criteria. The 90% C.L. upper limits for a downgoing flux of fast (f3 > 0.03) IMM's 
coming from above are at the level of 1.3 • 10~ 15 cm~ 2 sr _1 s _1 . The complete /3-dependence for MMs 
with g = go, 2gjj, 3gn and for dyons is shown in Fig. [71 




(1) 
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Figure 7: 90% C.L. upper limits for a downgoing flux of IMMs with g = go, *2gD, 3<?d and for dyons 
(M+p , g — go) plotted vs f3 (for strong etching). The poor limits at j3 ~ 10 -3 arise because the 
REL is below the threshold (for go and 2go) or slightly above the threshold (for 3gD and dyons), 
see Sect. 12.31 



Superheavy GUT magnetic monopoles in the cosmic radiation can traverse the Earth. Therefore the 
SLIM limit on their flux is one half of the IMM flux: 4> GUT < 6.5 • lCT 16 cm* 2 s" 1 sr" 1 for (3 > 0.03 for 
9 = 9d 0. 

Fig. [5] shows the flux upper limits for MMs of charge g = go and (3 > 0.05 vs monopole mass. 
Note that the SLIM limit is 1.3 • 10" 15 cm -2 sr* 1 s* 1 for MM masses smaller than ~ 5 • 10 13 GeV and 
0.65 • 10 -15 cm -2 sr _1 s _1 for masses larger than ~ 5 • 10 13 GeV. In Fig. [8] are also shown the limits 
obtained by the MACRO [B] and OHYA [53] experiments for g = go magnetic monopoles with j3 > 0.05. 

SLIM is the first experiment to extend the cosmic radiation search for Magnetic Monopoles to masses 
lower than the GUT scale with a high sensitivity. 

The addition of SLIM data to the MACRO data would improve the MACRO limits by only 18%. 

Large scale underwater and under ice neutrino telescopes (Amanda, IceCube, ANTARES, NEMO) 
have the possibility to search for fast IMMs with (3 > 0.5 to a level lower than the Parker bound [HI [25]. 
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Figure 8: Flux upper limits for cosmic MMs of charge g = go and (3 > 0.05 vs monopole mass. The 
figure shows the 90% CL. limits obtained by the SLIM, MACRO [6] and OHYA [23] experiments. 
MMs with masses smaller than ~ 5 ■ I0 13 GeV are detected only if coming from above; MMs with 
masses larger than ~ 5- 10 13 GeV can traverse the Earth, so an isotropic flux is expected. The Parker 
bound [26], obtained from the survival of the galactic magnetic field, and the limit obtained from the 
mass density for a uniform density of monopoles in the Universe |27| are also plotted. 
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